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Abstract 


impedance, kg/m'^ sec 


A time- dependent numerical formulation is de- 
rived for sound propagation in a two-dimensional 
straight soft- walled duct in the absence of mean 
flow* The time-dependent governing acoustic-dif- 
ference equations and boundary conditions are de- 
veloped along with the maximum stable time incre- 
ment. Example calculations are presented for sound 
attenuation in hard- and soft-wall ducts. The time- 
dependent analysis has been found to be superior to 
the conventional steady numerical analysis because 
of much shorter solution times and the elimination 
of matrix storage requirements. 
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cell coefficient 
cell coefficient 
cell coefficient 

ambient speed of sound, m/s 
cell coefficient 
cell coefficient 

frequency, Hz 
cell coefficient 
cell coefficient 

height of duct, m 

number of axial grid points 

number of transverse grid points 

length of duct, m 
transverse mode number 

time- dependent acoustic pressure, P*/p^c^^ 

o o 

spatially dependent acoustic pressure 

spatially dependent solution of Helmholtz 
equation 

period, 1/f", sec 
dimensionless time, t /T 
time step 

axial acoustic velocity, u /c 

o 

transverse acoustic velocity, ^ 

axial coordinate, x*/H^ 
axial grid spacing 

dimensionless transverse coordinate, y /h 
transverse grid spacing 
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e exrt condition 
i axial index (fig* 1) 
j transverse index (fig. 1) 
m cell index 
o ambient condition 
Superscripts : 

* dimensional quantity 
k time step 

(1) real part 

(2) imaginary part 


I, Introduction 


Both finite-difference and finite- element nu- 
merical techniques (refs. 1 to 27) have been devel- 
oped to study sound propagation with axial variations 
in Mach number, wall impedance, and duct geometry as 
might be encountered in a typical turbojet engine. 
Generally, the numerical solutions have been limited 
to low-frequency sound and short ducts, because many 
grid points or elements were required to resolve the 
axial wavelength of the sound. As shown in refer- 
ence 1 (eq, 77) for plane wave propagation, the num- 
ber of grid points or elements is proportional to 
the sound frequency and duct length, ^and inversely 
proportional to one minus the Mach number (ref, 2, 
fig. 6). This later dependence severely limits the 
application of numerical techniques for high Mach 
number inlets. 

Customarily, the pressure and acoustic veloci- 
ties are assumed to be simple harmonic functions of 
time,* thus, the governing linearized gas -dynamic 
equations (ref. 28, p. 5) become independent of time . 
The matrices associated with the numerical solution 
to the time independent equations must be solved ex- 
actly using such methods as Gauss elimination. It- 
eration techniques are unstable* As a result, large 
arrays of matrix elements must be stored idiich tax 
the storage capacity of even the largest computer. 

In unpublished work at Lewis using reference 29 as 
well as the work of Quinn (ref. 23, p. 3) the matrix 
has been modified to allow iteration techniques; un- 
fortunately the convergence is too slow to be of any 
practical value. Other approaches, such as in ref- 
erence 30, might still offer iterative possibilities . 



( 1 ) 


Some special techniques have been developed to 
overcome the above mentioned difficulties . As shown 
in references 3 and 10, the wave envelope numerical 
technique can reduce the required number of grid 
points by an order of magnitude* In reference 20, 
this technique was used to optimize multielement 
liners of long lengths at high frequencies. At the 
present time, this technique has been applied only 
to the simple cases of no flow and plug flow. A 
numerical spatial marching technique was also de- 
veloped in references 15 and 18. Compared to the 
standard finite- difference or finite- element bound- 
ary value approaches, the numerical marching tech- 
nique is orders of magnitude shorter in computa- 
tional time and required computer storage. The 
marching technique is limited to high frequencies 
and when reflections are small. 
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These and other symbols are defined in the nomen- 
clature. The dimensionless frequency r\ is de- 
fined as 
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As an alternative to the previously developed 
steady state theories, a time-dependent numerical 
technique is developed herein for noise propagation 
in a two-dimensional soft-wall duct in the absence 
of mean flow. Advantageously, matrix storage re- 
quirements are significantly reduced in the time- 
dependent analysis. The analysis begins with a 
noise source radiating into an initially quiescent 
duct. This explicit method calculates stepwise in 
real time to obtain the transient as well as the 
^steady” state solution of the acoustic field. The 
total time required for the analysis to calculate 
t|ie "steady" state acoustic field will determine the 
usefulness of the time -dependent technique. 

Time- dependent numerical techniques have been 
applied extensively to both one- dimensional sound 
propagation (ref. 31, p. 258), two-dimensional vi- 
bration problems (ref. 32, p. 452), and the more 
general problem of compressible fluid flow (ref. 33). 
References 31, 34, and 35 discuss in detail the sta- 
bility of numerical solutions to the wave equation. 
Herein, these teedmiques will be extended to include 
soft-wall impedance boundary conditions which would 
be encountered in inlets and exhaust ducts of turbo- 
fan engines . 

In the present paper, the governing acoustic- 
difference equation and the appropriate boundary 
conditions associated with time -dependent propaga- 
tion are presented. Next, the von Neumann method 
is used to develop the relationship between sound 
frequency and grid spacing to determine the maximum 
stable time increment. Immediately following the 
mathematical development, numerical solutions are 
presented for one- and two-dimensional hard- and 
soft-wall ducts. The results are compared with the 
corresponding steady analytical results. Finally, 
the time required to perform both the time-dependent 
and steady analyses are compared for increasing num- 
ber of grid points. 


Governing Squations and Boundary Conditions 

The propagation of sound in a two-dimensional 
rectangular duct, as shown in figure 1, is described 
by the linearised continuity and raomentuoi equations 
and the appropriate impedance boundary conditions. 


The asterisks denote dimensional quantities. 

The foregoing dimensionless equations apply to 
the scaled Cartesian coordinate system in which the 
height ranges between 0 and 1 and the dimensionless 
length ranges between 0 and L /h . 

Wave Equation 


Equations (1 to 3) now are combined to yield 
the dimensionless wave equation 




2 ^ 


ay2 


(5) 


Equation (5) in difference form will be solved to 
determine the pressure in the duct. 


Wall Boundary Condition 


The boundary condition at the surface of a 
sound absorbent soft-wall duct can be expressed in 
terms of a specific acoustic impedance defined as 





( 6 ) 


Substituting equation (6) into equation (3) yields 


By ^ Bt ^2 Bt 


(7) 


In the example problems to be considered, the impe- 
dance will be assumed constant with time. There- 
fore, equation (7) reduces to 


Bp 3 Bp 

a?= £ H 
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At the lower wall, the sign on g is changed to ac- 
count for the vector nature of v. It is also con- 
venient to express the specific acoustic impedance 
in terms of resistance 6 and reactance X as 


^ = 0 + iX 


(9) 


Entrance Condition 


Continuity and Momentum 

The lineariJted equations for mass and moraentum 
conservation can be written (ref. 28, p. 5) for a 
Cartesian coordinate system in the following dimen- 
sionless form: 


The boundary condition at the source plane 
P(o,y,t) can be of any general form with both trans- 
verse variations in pressure and multiple frequency 
content. However, the numerical technique will be 
compared later to previous solutions in which the 
pressure and acoustic velocities were assumed to be 
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plane waves at the entrance and to vary as 

or in dimensionaless form as l^ierefore 

the source boundary condition used here in 


where i and j denote the space indices 3 k the 
time index^ and Ay 3 and At are the space and 

time mesh spacing, respectively. All spacings are 
assumed constant. The time is defined as 


P(o,y,t) = ( 10 ) 


Exit Impedance 

In a manner similar to the wall impedance, the 
axial impedance at the duct exit can be defined as 


® u(L*/H*,y,C) 


( 11 ) 


= t + At =- (k + i)At 


(14) 


Solving equation .(13) for the pressure yields 

J 



For the plane wave propagation to be considered 
herein, is taken as 1 , which is exact for plane 
wave propagation in an infinite hard-wall duct. Al- 
so choosing to be 1 has lead to close agree- 

ment between numerical and analytical results for 
plane wave propagation into a soft-wall duct (refs. 
1 and 3). More general values for the exit impe- 
dance can be found in references 7 , 15 (eq. 

16, or 18 (fig. 7 ). * 

Initial Condition 

■ For times equal to or less than zero, the duct is 
assumed quiescent, that is, the acoustic pressure and 
velocrtres are taken to be zero. For times great— 
bhan zero, the application of the noise source 
(eq. ( 10 )) will drive the pressures in the duct. 



where cl is defined as 



C15) 


(16) 


Complex Notation 

Because of the introduction of complex no tar 
tion for the noise source and wall impedances, all 
the dependent variables are complex. The super- 
script ( 1 ) will represent the real term while ( 2 ) 
will represent the imaginary term; 

p = + ip(2) ^2) 

A siniilar notation applies to the acoustic velocities . 


15 an aigoritnm wnicti permits 
inarching out solutions from known values of pres- 
sures at times associated with k and k - 1. The 
procedure is explicit since all the past values of 
F are known as the new values of k -f 1 are com- 
puted. For the special case at t ^ 0, the values 
of the pressure associated with the k - 1 value 
are zero from the assumed initial condition. The 
parameter a was introduced into equation (15) be- 
cause it will play an important role in determining 
the stability (error growth) of equation (15) in 
this explicit iteration scheme. 


Difference Equations 

Instead of a continuous solution for pressure 
in space and time, the finite-difference approxi- 
mations will determine the pressure at isolated grid 
points in space as shown in figure 1 and at discrete 
time steps At. Starting from the known initial 
conditions at t = 0 and the boundary conditions, 
the finite-difference algorithm will march-out the 
solution to later times. 


Central Region fCell 


Away from the duct boundaries, in cell of 
figure 1 , the second derivatives in the wave equa- 
tion (eq, (5)) can be represented by the usual cen- 
tral differences in time and space (ref, 34 , p. 99 ) 
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^1=2 J 
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Ay^ 



( 13 ) 


Boundary Condition (Cells §2 to # 6 ) 

The expressions for the difference equations at 
the wail boundary are complicated by the impedance 
condition and the change in geometry of cells 
through #6 in figure 1. The governing difference 
equations can be developed by an integration process 
in which the wave equation (eq. ( 5 )) is integrated 
over the area of the cells and time: 


J ^t+At/2 „ 

f ff 


dx dy dt = 0 


(17) 

The spatial integration over the cell area is fully 
documented in reference 10, appendix D. In equa- 
tion (17) of this paper, the spatial integration 
method of reference 10 has been extended to include 
time. An illustration of the application of equa- 
tion (17) to cell ^2 is given in the appendix of 
this paper. 
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The finite-difference approximation for the 
various cells shown in figure 1 are expressed in 
terms of the cell coefficients a through g : 



method is often used to study the stability of the 
difference approximations to the wave equation. Ap- 
plication of the von Neuman method fref. 34, p* 104) 
to equation (13) requires that a in equation (16) 
be less than 1; which limits the time increment 



The derivation in reference 34 was for only one space 
dimension; however, the extension of the von Neuman 
method to two space dimension is relatively easy* 
When the time step satisfies equation (22), a condi- 
tion of linear instability exists (ref, 34, p. 106), 
which guarantees that a propagating acoustic mode 
can travel undirainished in a hard-wall duct (see 
discussion, ref. 18, p. 302). 



as) 


The subscript m denotes the cell number. These 
coefficients are listed in table I, 

Spatial Mesh Size 


Steady State Pressures 

In the sample problems to be presented in the 
next section, the time- dependent results will be 
compared to the results of the steady harmonic solu- 
tions of reference 10. The purpose of the section 
Is to show the rationale for constructing a steady 
state solution from the time-dependent results. 

Steady Harmonic Solution 

The steady harmonic pressure p (x,y) is de- 
fined as a solution to equation (5) when the pres- 
sure is assumed to be a simple harmonic function of 
time ; 


The mesh spacings Ax and Ay must be re- 
stricted to small values to reduce the truncation 
error. To resolve the oscillatory nature of the 
pressure the required number of grid points in the 
axial direction suggested in reference 1 was 


= Pg (23) 

Substituting equation (23) into equation (5), the 
wave equation takes the form of the classic Helm- 
holtz equation 


I 



(19) 


No requirement for the size of transverse spacing 
Ay was given in reference 1 other than the number 
of transverse grid points be increased until con- 
vergence is achieved. 

In the rectangular duct shown in figure 1, 
propagating transverse acoustic pressure modes (cos 
njry) can exist in the duct (ref. 15, eq. (B2)) when 

^ 2x\ (20) 

To resolve all the propagating modes, the number of 
grid points in the transverse direction, J, sugges- 
ted here is 


J 12ti (21) 

Equation (19) in conjunction with equation (21) 
would lead to equal axial and transverse mesh spac- 
ing which generally minimizes the truncation error 
(ref. 33, p. 288). 

Stability 

In the explicit time marching approach used 
here, round-off errors can grow in an unbounded 
fashion and destroy the solution if the time incre- 
ment At is taken too large. The von Neuman 


^ ^s ^ Ps 2 

2 + 2 P, = 0 (24) 

9x Sy 

In this case, where the source is a simple harmonic 
function of time, p represents the Fourier trans- 
form of P(x,y, t) (ref. 28, p. 11). The boundary 
conditions can also be modified by equation (23) as 
shown in reference 10. 

For a semiinfinite duct (or an equivalent fi- 
nite duct with P*c* exit impedance) with plane 
wave propagation, and hard walls, the solution for 
Ps is 


_ -i27Ttlx 


(25) 


In the next section, a transient solution to this 
problem will be compared to equation (25) . 


Transient Solution 


Recall, at the start of the numerical calcula- 
tion, the acoustic pressures and velocities were 
assumed zero throughout the duct and a pressure 
source begins a harmonic oscillation at x = 0 for 
t 2: 0. For the special case of plane wave propaga- 
tion in a hard-wall, semiinfinite duct, the analyti- 
cal solution to the wave equation (eq. (5)) is 
(ref. 36, p. 305) 
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P(x,t) 


^“i27mx gi^TTt 


0 < t < t)x (26) 
t > Tix (27) 


The pressure p(x>y) is now defined by divid- 
ding the instantaneous pressure P(x,y,t) by 
to obtain 


wall duct for plane wave propagation with ^ =1 

exit impedance (equivalent to semiinfinite duct) • 

The calculation was made with a length to height ra- 
tio (L*/H*) of 1 and a dimensionless frequency q 
of 1. The analytical and mimerical values of the 
acoustic pressure profiles along the duct are shown 
in figure 2, As seen in figure 2, agreement between 
analytical and the numerical theory is good. 


P(x,y) 


gi2irt 


<28) 


Using this definition, equations (26) and 
come 


fo 


P(x> 


-i27TT)X 


0 < t < qx 
t > qx 


(27) be- 

(29) 

(30) 


Consequently, the transient solution (eq* (30)) for 
the steady harmonic pressure equals the Fourier 
transform solution (eq. (25)) when 


Soft-Wall Duct 


As another example of the time- dependent analy- 
sis, the pressure distributions are computed for the 
case of plane wave propagation with a ^ =1 exit 

Impedance and a wall with impedance values of 0»16- 
i0*34. The calculation was made with a length- to- 
height ratio of 0,5 and a dimensionless frequency of 
0.6, The results of the time- dependent analysis 
along with the results of the solution of the equi- 
valent steady state Helmhomtz equation are d!is played 
in figure 3, The numerical results for the steady 
spatial solution p (x,y) are tabulated in appendix F 
of reference 10, As seen in figure 3, again the 
steady and time-dependent analyses are in good 
agreement* 


t > qx 


(31) Grid Point Variations 


In terms of real variables, equation (31) can 
be written as 


* 

* 

c 

o 


(32) 


The transient time t represents the time for the 
wave to travel down to the end of the duct, x*^ = L** 
Therefore, for the special case of one -dimensional 
plane wave propagation, the initial transient will 
pass when equation (31) holds. 


Since it may be desirable to integrate the wave 
with time to obtain a rms (root mean squared) pres- 
sure, the transient calculations will be continued 
into the steady domain for one period of oscilla- 
tion before the Fourier pressure p is calculated. 
Therefore, in this paper. 


and 


t = 
c 


= 31^+1 
H* 


p{x>y) 


_ P(x>y.l:) 


(33) 


(34) 


Figure 4 shows the effect of increasing the 
number of grid points on the computational time of 
the time- dependent approach for the hard -wall duct 
associated with figure 2, Roughly, as seen in fig- 
ure 4, the computational time is proportional to the 
number of grid points used. This is a considerable 
advantage over the steady technique in which the 
computational time more nearly increases with the 
square of the total grid points . 


Figure 5 shows a comparison of the computational ' 
times for the steady and time -dependent problems 
associated with the hard-wall duct shown in figure 2. 
As seen in this figure, the time- dependent analysis 
is considerably faster than the steady analysis. 

The value of J (grid points in y- direction) was re- 
stricted to 20 because of practical limitations on 
the size of the matrix which could be effectively 
handled in the steady analysis. When J was in- 
creased to 50 with 50 axial grid points, the steady 
analysis required 5500 seconds as compared to less 
than 20 seconds for the time-dependent analysis . 

This large increase in the steady state solution 
time results because of the manner In which the gen- 
eral matrix was partitioned (ref, 10, p. 14). In 
this case, the storage and computational times are 
proportional to the total number of transverse grid 
points squared (J)^. 


For more complicated problems, such as with higher- 
order acoustic pressure modes or where reflections 
are important, should be increased in succes- 
sive steps to check for convergence. 


Finally, these calculations were performed on 
the Univac 1100 computer. Faster computers can be 
expected to significantly reduce these calculational 
times . 


Sample Calculations 

In two sample problems to follow, the time- 
dependent results will be compared to the results 
of the steady harmonic solutions of reference 10. 

Hard-Wall Duct 

Numerical and analytical values of the pres- 
sure p(x,y) are computed for the case of a hard- 


Concluslons 

With the possible exception of the wave enve- 
lope technique (ref. 10) or the spatial marching 
technique (ref, 18), the numerical time- dependent 
method of analysis represents a significant advance 
over previous steady numerical theories. By elimi- 
nating large matrix storage requirements, numerical 
calculations of high sound frequencies are now pos- 
sible. Also, because matrix manipulation is not re- 
quired, the time- dependent approach is simpler to 
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program and debug. Although flow has not been con- 
sidered herein, the extension to the more general 
flow situation appears to be 'straightforward. 


APPENDIX 

Finite-Difference Equations and Coefficients 


where 



The derivation of the difference equation for 
the various cells starts with equation (17), which 
is rewritten here as 


cm 


2 ^ 




Sy- 


^Idx dy dt = 0 


(Al) 


where the plus sign (+) in the upper limit of inte- 
gration means to evaluate the parameters along 
either the upper or right-hand boundary of the in- 
tegration cell, shown in figure 1 by the dashed 
lines, while the negative sign (-) applies to either 
the lower or left-hand boundary of the integration 
cell, depending on whether x or y is consid- 
ered. The (+) and (-) notation was used since the 
spatial integration limits will vary from cell to 
cell. For cell for example, in the x inte- 
tral, (+) is represented by x + Z^/2 and (-) by 
X - Ax/2, while in the y integral, (+) is repre- 
sented by 1 and (-) by 1 - Ay/2* 

The pressure P can be assumed constant over 
the cell area and likewise the second derivative in 
time can also be assumed constant over the entire 
cell area. Therefore, moving the spatiai deriva- 
tives to the right side of equation (AI) gives 


/*t+At/2 o 

^ 

*/t-At/2 




dt 


pt^At/2 

. / 

t-At/2 


X 



(A2) 



(A5) 


The parameter can now be evaluated for each 

cell (labeled ra) in figure 1. The procedure for 
each cell is given in detail in reference 10, appen- 
dix D. To illustrate how the time integration Is 
coupled to the space integration, Che derivation for 
cell #2 will now be presented. 


Por cell #2, 



dy dx = ^ 


(A6) 



P*' 

_ 1+1. i 


- 2P^ . + P^ . 

iiJ k-_1.3 




(A7) 


i+ 


Sp 

dy 


dy 


£ dt 


' P" - P^ . , , 

JjJ iUzi 

Ay 


CA8) 


The expression for 


in equation. (A8) is the 


wall boundary condition as given by equation (8) in 
the body of this report. Substituting equations 
(A6) , (A7), and (A8) into equation (A5) yields 


but 


L 


t+At/2 


t-At/2 




t+At/2 / ^1 . pk 

= l_iiJ iil 

I At 
t-At/2 \ 


(pk _ pk-l\ pk+1 . 2P>^ . + pk-1 

1.3 J->J = _iiJ liJ LiJ. 

At i At 


Therefore, equation (A2) can be written as 


2 Ay 


2 At Ax 


- 2P^ + p^ 

_i+LJ ** ^ ^ 


t+At/2 


i-1.1 f 

^ ^*4-At/2 


t+At/2 


(A3) 


X dt 


X dx 


r+ „ r" "'"hs. . r+ 

f ^ J dt j 

t-At/2 


/ k k \ .+ 

('l.J - "‘I 


(A9) 





(A4) 


Equation (A9) can be further simplified by 
noting 
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t+At/2 


t-At/2 


dt = At 
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7. Sigman, R. K. , Majjigi, R, K. , and Zinn, B. T. , 
"determination of Turtofan Inlet Acoustics 
Using Finite Elements," AlAA Journal, Vol. 16. 
Nov. 1978, pp. 1139-1145. 
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t+At/2 


t-At/2 


l|dt = P. . 

Ot 1, J 


t+At/2 

t-At/2 






CA13) 


Substituting equations (AlO) to (Al3) into equa- 
tion (A9) yields 




2 



+ 2P 


k 

j-1 



+ 



P 


k 

i+l^j 


. a ^ k+i ji ^ k-i 
5 At ^ 5 At ^i,j 


(A14) 


which contain the coefficients that appear in 
table I for cell #2. The coefficients for the other 
cells are found in a similar manner. 


8. Abrahamson, A. L. , "A Finite Element Algorithm 

for Sound Propagation in Axisyrnmetric ducts 
Containing Compressible Mean Flow," Wyle 
Labs. , Inc. , Hampton, VA, June 1977. (NASA 
CR-145209) 

9. Craggs, A. , "A Finite Method for Modelling Dis- 

dipative Mufflers with a locally Reactive 
Linning," Journal of Sound and Vibration . 

Vol. 54, Sep. 1977, pp. 285-296. 

10. Baumeister, K. J., "Finite-Difference Theory for 

Sound Propagation in a Lined Duct with Uni- 
form Flow Using the Wave Envelope Concept," 
NASA TP-1001, 1977, 

11. Kagawa, Y. , Yamabuchi, T. , and Mori, A. , "Finite 

Element Simulation of an Axisymma trie Acoustic 
Transmission System with a Sound Absorbing 
Wall," Journal of Sound and Vibration. Vol. 
53, Aug- 1977, pp- 357-374. 

12. Eversman, W. , Astley, R. J. , and Thanh, V. P. , 

"Transmission in Nonuniform Ducts - A Compar- 
ative Evaluation of Finite Element and 
Weighted Residuals Computational Schemes," 
AIAA Paper 77-1299, Oct. 1977- 


References 

1. Baumeister, K. J. and Bittner, E. C. , "Numerical 

Simulation of Noise Propagation in Jet Engine 
Ducts," NASA TN D-7339, 1973. 

2. Baumeister, K. J. and Rice, E. J. , "A Difference 

Theory for Noise Propagation in an Acoustically 
Lined Duct with Mean Flow, " Aeroacoustics ; Jet 
and Combustion Noise; Duct Acoustics, Progress 
in Astronautics and Aeronautics, Vol. 37, Aar- 
erican Institute of Aeronautics and Astronau- 
tics, New York, 1975, pp. 435-453. 

3. Baumeister, K. J. , "Analysis of Sound Propaga- 

tion in Ducts Using the Wave Envelope Con- 
' cept," NASA TN D-7719, 1974. 

4. Quinn, D. W. , "A Finite Difference Method for 

Computing Sound Propagation in Nonuniform 
Ducts," AIAA Paper 75-130, Jan. 1975. 

5. Baumeister, K. J. , "Wave Envelope Analysis of 

Sound Propagation in Ducts with Variable Axial 
Empedance," Aeroacoustics: Duct Acoustics, 

Fan Noise and Control Rotor Noise, I. R. 
Schwartz, H. T. Nagamatsu, andW. Strahle , eds.'. 
Progress In Astronautics and Aeronautics, Vol. 
44, American Institute of Aeronautics and As- 
tronautics, New York, 1976, pp. 451-474. 

6- Quinn, D. W, , "Attenuation of Sound Associated 
with a Plane Wave in a Multisectional Duct," 


13. Watson, W. R. , "A Finite Element Analysis of 

Sound Propagation in a Rectang ular Duct of 
Finite Length with Peripherally Variable 
Liners," AIAA Paper 77-1300, Oct. 1977. 

14. Abrahamson, A. L., "A Finite Element Algorithm 

for Sound Propagation in Axisymmetric Ducts 
Containing Compressible Mean Flow, " AIAA Paper 

77- 1301, Oct. 1977. 

15. Baumeister, K. J. , "Numerical Spatial Marching 

Techniques for Es tintat in g Duct Attenuation and 
Source Pressure Profiles," 95th Meeting 
Acoustical Society of America, Providence, 
Rhode Island, May 16-19, 1978- (Also NASA TM- 
78857, 1978) 

16. Tag, I. A. and Lumsdaine, E. , "An Efficient 

Finite Element Technique for Sound Propagation , 
in Axisymmetric Hard Wall Ducts Carrying High 
Subsonic Mach Number Plows, AIAA Paper 

78- 1154, July 1978. 

17. Astley, R. J- and Eversman, W. , "A Finite El^ 

ment Method for Transmission in Non-Uniform 
Ducts without Flow: Comparison with the 

Method of Weighted Residuals," Journal of 
Sound and Vibration, Vol- 57, Apr. 1978, 
pp. 367-388. 

18. Baumeister, K. J. , "Numerical Spatial Marching 

Techniques in Duct Acoustics," Journal of the 
Acoustical Society of America, Vol. 65, Peb. 
1979, pp. 297-306. 


7 



19. Majjigi, R. K. , "Aj>pllcation of Finite Element 

Techniques in Predicting the Acoustic Proper- 
ties of Turhofan Inlets," Ph.D. Thesis, 

Georgia Institute of Technology, Atlanta, GA, 
1979. 

20. Baumeister, K. J. , "Otpimized Multisec tioned 

Acoustic Liners," AIAA Paper 79-0182, Jan. 

1979. 

21. Majjigi, R. K. , Sigman, R. K. > and Zinn, B. T. , 

"Wave Propagation in Ducts Using the Finite 
Element Method," AIAA Paper 79-0659, Mar. 1979. 

22. Astley, R. J. and Eversman, W. , "The Application 

of Finite Element Techniques to Acoustic Trans- 
mission in Lined Ducts with Flow," AIAA Paper 
79-0660, Mar. 1979. 

23. Quinn, D. W. , "A Finite Element Method for Com- 

puting Sound Propagation in Ducts Containing 
Flow," AIAA Paper 79-0661, Mar. 1979. 

24. Abrahamson, A. L. , "Acoustic Duct Liner Optimi- 

zation Using Finite Elements," AIAA Paper 
79-0662, Mar. 1979. 

25. Tag, I. A. and Akin, J. E. , "Finite Element So- 

lution of Sound Propagation in a Variable Area 
Duct," AIAA Paper 79-0663, Mar. 1979. 

26. Lester, H. C. and Parrott, T. L. , "Application 

of Finite Element Method for Computing Grazing 
Incidence Wave Structure in an Impedance Tube; 
Comparison with Experiment," AIAA Paper 
79-0664, Mar. 1979. 

27. Baumeister, K. J. and Majjigi, R. K- , "Applica- 

tions of Velocity Potential Function to Acous- 
tic Duct Propagation and Radiation from Inlets 


Using Finite-Element Theory," AIAA Paper 
79-0680, Mar. 1979. 

28. Goldstein, M. E. , Aeroacoustics , McGraw-Hill, 

New York, 1976. 

29. Beaubien, M. J, and Wesler, A, , "Iterative, 

Finite Difference Solution of Interior Eigen- 
values and Eigenfunctions of Laplace, s Opera- 
tor," Computer Journal , Vol. 14, Aug. 1971, 
pp. 263-269. 

30. Browne, B. T. and Lawrenson, P. J. , Numerical 

Solution of an Elliptic Boundary- Value Problem 
in the Complex Variable," Institute of Mathe- 
matics and its Applications Journal , Vol. 17, 
1976, pp. 311-327. 

31. Richtmyer, R. D. and Morton, K. W. , Difference 

Methods for Initial-Value Problems , 2nd ed. , 
Inters cieace, New York, 1967. 

32. Gerald, C. F. , Applied Numerical Analysis , 2nd 

ed. , Addis on- Wesley, Reading, MA, 1978. 

33. Roache, P. J. , Computational Fluid Dynamics , 

Hermosa, Albuquerque, NM, 1972. 

34. Clark, M. and Hansen, K. F. , Numerical Methods 

of Reactor Analysis , Academic Press, New York, 
1964. 

35. Hildebrand, P, B. , Methods of Applied Mathema- 

tics , Printice-Hall, Englewood Cliffs, NJ, 
1952. 

36. Budak, B. M. , Samarskii, A. A, , and Tikhonov, 

A. N. , A Collection of Problems on Mathemati- 
cal Physics , Pergamon, (Oxford), 1964. 


TABLE I. - COEFFICIENTS IN DIFFERENCE EQUATIONS 



8 





Figure 1. - Grid-point representation of two-dimensional flow duct. 
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Figure 2 . - Anaiytical and numerical pressure profiles for 
one-dimensional plane wave sound propagation in a hard 
wall duct for 17 = 1 and L*/H* = 1 , 
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Figure 3. - Pressure profiles on the duct axis for incident 
plane wave sound propagation in a soft wall duct (?? = 0 6 
L-/H‘’-==0.5,?^^ = 0.16-10.34). 
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Figure 4. - Effect of increasing number of grid points on 
calculational time of transient solution for plane wave 
propagation in a hard wall duct ( 7 ? = 1, L*/H* = 1). 
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Figure 5. - Comparison of transient and steady state calcu- 
lations for plane wave propagation in a hard wall duct 
(T? = 1, L*/H* = 1, J = 20). 
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